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Abstract. We develop a theoretical model to study the Intensity-Intensity 
correlation of polarisation entangled photons emitted in a biexciton-exciton cascade. 
We calculate the degree of correlation and show how polarisation correlations are 
affected by the presence of dephasing and energy level splitting of the excitonic 
states. Our theoretical calculations arc in agreement with the recent observation of 
polarisation dependent Intensity-Intensity correlations from a single semiconductor 
quantum dot [R. M. Stevenson et. at Nature 439, 179 (2006)] . Our model can 
be extended to study polarisation entangled photon emission in coupled quantum dot 
systems. 
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1. Introduction 

Polarisation correlations of photons emitted in cascade emission are well known 
phenomenon and numerous theoretical and experimental studies exist in the literature 
on this subject since the early days of quantum optics [1, 2, 3, 4, 5, 6]. Some of the 
earlier studies were motivated in testing generalised Bell's inequalities [2] , the existence 
of hidden variables and whether quantum mechanics was a non-local theory or not [3, 4], 
following the question raised by Einstein, Podolsky and Rosen [7]. In recent times 
polarisation correlated photon pairs have become important in the field of quantum 
information science due to their entangled nature. Moreover many applications of 
quantum information, such as quantum key distribution [8], efficient optical quantum 
computing [9], long distance quantum communication using quantum repeaters [10] and 
implementation of quantum telecommunication schemes [11] require single photon pairs 
per cycle. This requirement of entangled photon pairs per cycle of excitation could be 
satisfied by cascade emission from a single atom or atom like systems like semiconductor 
quantum dots, provided one gets over the inherent asymmetries. Recently such cascade 
emission has been reported for semiconductor quantum dots [12, 13, 14, 15]. It was 
further seen that polarisation entanglement of the emitted photon pairs was degraded 
by the presence of energy level splitting of the intermediate cxcitonic states and any 
incoherent process that leads to a population transfer between the two intermediate 
excitonic states [14, 15, 16]. Moreover dephasing arising due to interaction of the 
quantum dot with its solid state environment can also degrade the entanglement 
[17]. Some recent studies have also shown how the fidelity of entanglement depends 
on excitonic level splitting [18] and the dynamics of the incoherent dephasing [19]. 
Different methods have been proposed to reduce and control the incoherent dephasing 
and energy level splitting of the excitonic states thereby preserving the entanglement 
in the system [14, 15, 16, 18, 20, 21, 22]. Further, methods to enhance the generated 
entanglement by coupling the quantum dot to a micro-cavity have also been proposed 
[23, 24]. As quantum dot systems arc of great importance for future applications in 
quantum information science, a clear yet simple model for understanding the effects of 
all these different decoherence mechanism on the dynamics of the system is required. 
Thus we develop, in this paper, a simple theoretical model to analytically study the 
influence of different decoherence mechanisms and the intermediate state splitting on 
the generation of polarisation entangled photon pairs in cascade emission. 

2. Model 

We consider a four level system under going cascade emission as our model. We show 
a schematic diagram of such a cascade in flgure 1. The excited state \i) and the 
intermediate states \a), would correspond to the biexcitonic and optically active 
excitonic states respectively in a quantum dot. Further is taken to be the ground 
state. Here 27 = 2(71+73) is the total spontaneous emission rate of the state 272, 274 
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are the spontaneous emission rates of the states \a) and \[3) respectively and 27/3a(27a/3) 
is the incoherent dephasing rate of the state |a) The energy level splitting of the 

intermediate state is given by A. In this type of four-level cascade scheme there are two 




Figure 1. Schematic diagram of a four level cascade system. Here H and V 
refers to horizontally and vertically polarised photon emission. A is the energy level 
separation of the intermediate states and 7's are the spontaneous emission rates given 
by 7fc = 2ujf.i\dki\'^ /ihc^ . The incoherent dephasing rates of the intermediate states 
are given by 2^ap and 2'ypa respectively. 



decay paths for the excited state, \i) \a) — > |j) and \i) \P) The generation 

of entanglement in these scheme is attributed to the fact that this decay paths can 
become indistinguishable. The eigenbasis of this system is formed by the four states 
({K)}) {|<^)}) {|/5)}> {|i)})- this basis the radiative transition from the excited state 
generates collinearly polarised photons with linear polarisations along two orthogonal 
directions denoted by H (horizontal) and V (vertical). When the states |a) and \(3) are 
degenerate, the decay paths become indistinguishable and we get a maximally entangled 
two photon state [7, 12] 

\E) = -^ilHiH,) + IV^V^)). (1) 

In practical systems like atoms and quantum dots these levels are usually non degenerate 
and hence the entanglement of the emitted photon pairs depend completely on the degree 
of degeneracy and dynamics of these intermediate states. In our model we have taken 
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them to be non-degenerate and study the effect of such intermediate level splitting on 
the correlation of the emitted photon pairs. To understand the effect of incoherent 
dephasing and energy level splitting of the excitonic state on the dynamics of emitted 
photon pairs from the cascade, we need to study the two time second order correlations. 
This is given by, 

(//) = (e7,^,^^) ■ i-(f,t)e^,^,<^^) ■ E-{r,t + r) 

■ ^(02,02) ■ t + r)e(e^,^^) ■ E+(f, t)). 

(2) 

where {II) stands for the two time polarisation angle dependent intensity-intensity 
correlation (/(^j ,j(,2)(r, t + r)/(e^ ,^-^)(r, t)). Further E~^ {r, t) {E~ {r, t)) is the positive 
(negative) frequency part of the quantized electric field operator at a point r in the 
far-field zone and e(e,0) is the polarisation unit vector of the measured radiation at the 
detector along any arbitrary direction given by {0,(f)). e(e^(^)'s are related to the linear 
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Figure 2. Schematic diagram for orientation of the polarisation unit vectors. Here 
en and ey corresponds to horizontal and vertical polarisation respectively, e^^^ and 
e^^-* arc two arbitrary orthogonal pair of polarisation unit vectors. 



polarisation unit vectors in, ey (where H stands for horizontal and V for vertical) by 

cos^^ 



-e'"^ sin I 



e ^'^ sin 
cos 6* 



(3) 



and these satisfy the relation (e|g''^-) • ej^''^)) = 0. The above relation can be understood 
as an unitary transformation between a basis defined by the linear polarisation unit 
vectors and a basis defined by e^^^ and e*^^\ In experimental setup the angles ^, would 
correspond to the orientation of the optic axis of a half/quarter wave plate to the 
direction of propagation of the emitted radiation. Let us now consider for simplicity 
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that both the levels \a) and in figure. 1 have the same incoherent dephasing rates 
7a/3 = 7/3a- Further we assume that the spontaneous decay rates of the intermediate 
levels are also equal. Such assumptions are well justified as they do not infiuence the 
dynamics of the system significantly and yet leads to a simplified form of the second 
order correlation, thereby providing a better understanding of the problem. Under the 
above assumptions and for 0i = 02 = the form of the two-time polarisation angle 
dependent intensity-intensity correlation is found to be , 

{i{e2,t + T)m,t)) = ^,vlvj{\z){z\,) 

X + cos 201 cos 2^26-^^^^+^^"'^)^ 

+ sin 2^1 sin 2^26"^^^'+^"'^^^ cos( At)} (4) 

where Vi = \dai\ = \dpi\ and V2 = \dja\ = \dji3\. The above simple form of the 
second order correlation has been derived to match our theoretical analysis to that of 
the experiment results [15]. For details of the mathematical analysis leading to the 
generalised form of the two time intensity-intensity correlation the reader is referred to 
section 4 of this paper. One can clearly see from equation (4) that the second order 
correlation is profoundly infiuenced by both the incoherent dephasing rates as well as 
the energy level splitting of the intermediate states. Note further, that in the presence of 
small A this becomes equivalent to the second order correlations measured in ref.[3, 4]. 
Next we define a quantity the degree of correlation as, 

where /j,, jj,' stands for mutually orthogonal polarisation basis. The degree of correlation 
varies between -|-1 and —1, where +1 represent perfect correlation (—1 for anti- 
correlation) and represent no polarisation correlation. 

3. Results and Discussion 

3.1. Effect of excitonic level splitting on the correlation 

In figure. 3(a) we show how the time averaged degree of correlation varies with the 
basis angle for different values of splitting A, of the excitonic levels. Note that here 
the excitonic level dephasing 7^^ has been taken to be zero. We see that the degree 
of correlation is independent of the polarisation basis when A = and takes a value 
C/j, = 1. This correspond to perfect polarisation correlation among the emitted photons. 
From the expression of c^^ it is clear that this can happen only when the cross-polarised 
correlations vanishes and the emitted photons are perfectly co-polarised. One can even 
see this exphcitly from equation (4) by putting the values oi 61,62 — O1 + 7r/2 for H-V , 
D-D' and V-H basis, where H, V, D and D' stands for horizontal, vertical, diagonal and 
orthodiagonal polarisation basis respectively. Further as the cross-polarised correlations 
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Figure 3. (Colour online) (a)Degree of correlation averaged over time as a function of 
basis angle, for different excitonic level splitting A. H , D, D', V stands for horizontal, 
diagonal, orthodiagonal and vertical polarisation basis respectively. Here we have 
considered zero dephasing of the excitonic states, (b) Conditional measurement of 
intensity-intensity correlation in the circular basis. The red curve corresponds to co- 
polarised (6*1 — 62 = 7r/4, — 4'2 ^ — 7r/2)photons and the blue for cross-polarised 
(6*1 = 02 ~ 7r/4, 01 — — 7r/2,(/)2 — vr/2) ones. The solid curve is for A = and the 
broken one for A ~ large. Here R and L stands for right and left circular polarisation. 
The R-R correlation curve in case of large splitting is time shifted for better comparison 
to the R-L correlation. All parameters are normalized with respect to 7. 



are absent the pair of photons emitted in one excitation cycle can take either of the two 
paths \i) \a) \j) or \i) \/3) \j) thus making these paths indistinguishable. 
As a consequence we do not get the "Welcher Weg" or which path information thereby 
making the final state of the emitted photon pair entangled in both the linear and 
diagonal polarisation basis. The generated entangled states can hence be written as 
1/V2{\HH) + \VV)) and 1/V2{\DD) + \D'D')) for the rectihnear and diagonal basis 
respectively. Note further that in this case perfect anti-correlation (c^ = -1) is expected 
for measurement in the circularly polarised basis with the entangled state given by 
l/\/2(\RL) + \LR)). Thus one should get perfectly cross-polarised photons as the co- 
polarised correlations vanish in this basis. This is exactly what we get from the general 
expression of equation (4) [see section 4, equation (13)] and is shown by the solid curves 
in figure. 3(b). Further in figure. 3(a) we see that the degree of correlation is practically 
independent of the excitonic level splitting A in the rectilinear basis. As we change our 
polarisation basis the effect of A becomes significant. In the diagonal basis for example 
with the increase in level splitting the degree of correlation gradually decreases and 
eventually vanishes. In presence of A, the cross-polarisation does not vanish and we 
have a which path information for the emitted photons when we measure the second 
order correlations, thus destroying any entanglement in the system. The behaviour of 
the correlations in the circular basis in presence of large excitonic level splitting is shown 
by the broken curves in figure. 3(b). One can clearly see that there is no polarisation 
correlation at all for large A. The sinusoidal behaviour of for non zero value of A 
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as seen in figure. 3(a) is in agreement with the classical linear polarisation correlation 
behaviour. Note that our theoretical results are in agreement to experimentally observed 
data [15]. 

It should be noted that in our analysis we have concentrated on the calculation of the 
quantum correlation c^. This also was measured in the experiment of Stevenson et. 
al. We have not examined measures of entanglement like concurrence. This is because 
if A -the intermediate state exciton splitting is nonzero then horizontal and vertical 
photons have different frequencies-which amounts to saying that we have for nonzero 
A quantum states which are characterized by two different parameters and measures of 
entanglement in such situations do not exist. 



3.2. Effect of decoherence on the correlation 

In figures. 4 (a) and (b) we show how the incoherent dephasing of the intermediate 
excitonic states affect the time averaged degree of correlations when excitonic states 
are non-degenerate (A ^ 0) and degenerate (A = 0) respectively. Note that here we 
have assumed that both the intermediate states have same dephasing rates. One can 
clearly see that the affect is different for different measurement basis. The degree of 
polarisation correlation for example in the rectilinear basis decreases with increasing 
dephasing irrespective of whether the excitonic states are non-degenerate or degenerate. 
For large dephasing rates the emitted photon pairs become almost un-correlated in their 




Figure 4. (Colour online) (a)Degree of correlation averaged over time as a function of 
basis angle, for large A and different incoherent dephasing rates ^a[3 of the intermediate 
level. Here we have assumed that both the intermediate states dephase at same rate 
i-e 7a/3 = "f/3a- H , D, D', V stands for horizontal, diagonal, orthodiagonal and vertical 
basis respectively, (b) Same as (a) for A = 0. 



polarisation. This is attributed to the presence of significant cross-polarised correlation 
for large dephasing rates of the intermediate states. The incoherent dephasing of the 
intermediate levels causes a incoherent population transfer among the states \a) and 
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Figure 5. (Colour online) Conditional measurement of intensity-intensity correlation 
in the circular basis for incoherent dcphasing /l = 10. The red curve corresponds to 
co-polarised {9i = 62 = 7r/4, (pi ^ (f>2 ^ — 7r/2)photons and the blue for cross-polarised 
{61 = 02 ~ 7r/4, 01 = — 7r/2,02 — 7''/2) ones. The solid curve is for A ^ large and 
broken one for A = 0. Here R and L stands for right and left circular polarisation. The 
R-R correlation curve are time shifted for better comparison to the R-L correlation. 



1(3) thereby allowing the second photon to be emitted with orthogonal polarisation to 
the first one. In the diagonal basis on the other hand the dephasing does not affect the 
correlation at all for large A but significantly decreases the correlation when A = for 
large dephasing rates. So we see that in diagonal basis even when the intermediate levels 
are degenerate we can still have significant cross-correlation if there is some incoherent 
relaxation process by which they can get coupled. This in turn spoils the quantum 
correlation in the system as can be seen clearly from figure. 4(b). In figure. 5 we show 
how the correlations behave in the circular basis in presence of large dephasing rate 
(70/3/7 = 10) for both non-degenerate and degenerate intermediate states. We find that 
in the circular basis decoherence arising due to the incoherent dephasing does not affect 
the degree of correlation of the emitted photons when A is large. Further we find that 
for degenerate intermediate levels, even though the degree of correlation = — 1 for 
zero time delay, it vanishes at all later time in presence of the large dephasing. Thus 
the decoherence makes the perfectly anti-correlated photons completely uncorrelated. 
The incoherent relaxation process discussed by us here are practically present in the 
biexcitonic-excitonic cascade in quantum dots [17, 19]. Thus we have shown by a simple 
model how the decoherence arising due to such incoherent processes would strongly 
affect the quantum correlations. 

4. Detail derivation of the Intensity-Intensity correlation 

Our model consist of a biexcitonic state and two excitonic states labelled as \i) and 
\P) respectively. The equihbrium state is labelled as The biexcitonic state 
decays by emission of either a horizontally(II) polarised photon \a)) or a 

vertically(V) polarised photon — » |/9))[12, 13]. The excitonic state |q;)(|/?)) decays 
to the equilibrium state \j) by emission of a II(V)- polarised photon. The two excitonic 
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states have a energy difference of ^A. Note tliat tlie splitting of the excitonic state 
in quantum dots arises due to anisotropic electron- hole exchange interactions [25, 26]. 
Figure (1) show a schematic diagram of our model. The eigenbasis of this system 
is formed by the four states ({K)}, {!«)}, In this basis the unperturbed 

Hamiltonian 7i is given by, 

n^Y.^k\k){k\ (6) 

k 

Where ?iWk is the energy of the four levels {k — i,a, (3, j). Note that this kind of energy 
level scheme has been extensively used to study the cascade emission in quantum dots 
[14, 16, 17, 18, 19]. The spontaneous emission and dephasing effects in the system are 
incorporated via a master equation technique [27] under the Born, Markov and rotating 
wave approximations and is given by, 

£p = - -^{Sii, p} - 72{5'«a, P} - 1a{Si3!3, P} - ll3a{Saa, P} 
- 7a/3{'5/3/3,P} + 2(7lP + l3PiiSj3l3 + l2PaaSjj 

+ liP0pSjj) + 2 {'JpaPaaSpp + la/SPpisSaa) ■ (7) 

Here 27 = 2(71 + 73) is the total spontaneous emission rate of the biexcitonic state 
272( 274) and 27^30, (27q,;3) are the spontaneous emission rate and incoherent dephasing 
rate of the excitonic state |q;)(|/5)) (see Fig. 1). Such incoherent dephasing arises 
in quantum dots due to it's interaction with the solid-state environment (in form of 
spin flip processes or phonon scattering) [17]. The curly bracket {..,..} stands for 
the anti-commutator and Ski = \k){l\{Sli = \l){k\ ) is the atomic lowering (raising) 
operator which follow the simple angular momentum commutation relations. To study 
the dynamical evolution of this four-level cascade system we solve for the time evolution 
of the density operator which is given by, 

(8) 

We will set the energy of the state | j) equal to zero henceforth. On substituting equations 
(6) and (7) in (8) and solving for the population terms we get, 

Pa/3(i)=e-(«°-^^)Va/9(0); 

Paa{t) = e-''"* (^cosh(At) + ^sinh(At)) p„«(0) 

+ 2e-«°*^sinh(Ai)p^;3(0) 

+ ^C(t)+p,(0)e-2^*D(t); 
27 

Ppp(t) = e-"°* (cosh(At) - ^ smh{At)j pppiO) 
+ 2e-«°*^sinh(^t)pc«a(0) 

+ ^F{t)+p,{0)e-''''K{t). (9) 
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Here R signifies a constant feeding of population into the state \i) from some arbitrary 
state \n). Note that in our model we are only concerned with the dynamics of the 
cascade decay once the upper level is populated, and thereby do not consider explicitly 
the pumping of the biexciton state \i). Further the excitonic level splitting A 



, do = (72 + 74 + lap + 7/3a), = (74 - 72 + 7a/3 - 7/3a) and A = ^JtI + A'^ocHliSa- The 
time dependent coefficients C,D,F and K are given by, 



C{t) 



271(1 + ^) 



A), 

r 



(271 (1 + 
+ {A^- 

(273 (1 - 
+ {A^- 

(273 (1 - 
+ {A^-A). 



I , 737a/3 

AJ A 

A), 

, 7l7/3a 

AJ A 

A), 

AJ A 



\ _ g-(ao-A)t 

2(ao - A) 

I _ g-(ao-A-27)t 

2(ao- A-27) 

\ — g-(ao-A)t 

2(ao - A) 

\ — g-(ao-A-27)t 

2(ao - A - 27) 



(10) 



The effect of non-degenaracy of the excitonic states and their incoherent dephasing on 
the dynamical evolution of the system shows up if one studies the two-time nonclassical 
second order correlation defined in equation (2). For our four level system the explicit 
form of the positive frequency part of the electric field operator is given by [27] , 



+ 
+ 
+ 



fix {fix d,3i)\ \l3){i\t 

— * 1 

h X {fix dja) \j){a\t 
fi X {fi X djp) 



111 



\3)m- 

Finally using equation (11) in (2) we get the general form of the two time intensity- 
intensity correlation 

(//) = ^{([(e^ • 4)*cos^i|z)(a|, 

+ (eV-4)*e''^'sin^i|i)(/3|i] 

X {\eH ■ d cos^ 92\Q){a\t+r + ley ■ djp\^ sw? 92\(3){(3\t+r 

+ e~"^^{eH ■ d,ja)*{ev ■ dj^) cos sin 6'2|a) 

+ e"^^{eH ■ dja){ev ■ dj/s)* cos 62 sin $21(3) {a\t+r) 



X 



[in ■ dai) cos6i\a) {i\t + {iv ■ dpi) e sin6'i|/3)(i|t )}. 



(12) 
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The two time correlation function that appears in equation (12) is evaluated by invoking 
the quantum regression theorem [28] and equation(9). Finally we get, 



X {f,{T)+W,{T)+f2{T)+W2{T) 

+ (cos2^i + cos2^2)(/i(r) -u>2(r)) 

+ (cos 2^1 - cos2^2)(wi(r) - /2(r)) 

+ cos 2^1 cos 2^2 (/i(r) + u;2(r) - /2(r) - Wi{t)) 

+ sin2^isin2^2(e"*^*^+'^^^M(r) + e*(*i+^2)M*(r))}. 



(13) 



Here T>i — \daj\ — \dpj\ and T>2 — \dja\ — \djp\. The /'s, w's and u are found from the 
solutions of the density matrix equations (9) and are given by, 

/i(r) = 6""°^ (^cosh(Ar) + ^ sinh(Ar) 

/2(t) =2e-"°^^sinh(AT), 

w;i(r) = 2e-"°"^sinh(ylT), 

W2{r) = e-"°^ (cosh(AT) - -^sinh(^T) 
u{r) =e-("°-^^)^. (14) 

The equation (13) gives the most general form of the two time intensity-intensity 
correlations for arbitrary polarisation directions and for any system undergoing a cascade 
emission. In the special case of (pi — (f)2 — 0, 7q,/j = ^y^a ^-nd 72 = 74 this reduces to the 
simplified result (4) of section (2). 

5. Conclusions 

In conclusion we have developed a simple theory to understand how the dephasing 
and energy level splitting of the excitonic states can affect polarisation entanglement 
of photons emitted in a biexciton-exciton cascade. We have also shown how these 
effects are important in determining whether the emitted photon pairs are classically 
correlated or entangled in different polarisation basis. Further we have shown that our 
theoretical calculation is in agreement with the experimental results found in context to 
such cascade emissions in quantum dots. As a future prospect it would be interesting 
to extend the method of the present paper to a system of coupled dots or a dot in the 
micro-cavity. 
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